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Abstract 

Studies of limits on active-sterile neutrino mixing derived from big bang 
nucleosynthesis considerations are extended to consider the dependance of 
these constraints on the primordial deuterium abundance. This study is moti- 
vated by recent measurements of D/H in quasar absorption systems, which at 
present yield discordant results. Limits on active-sterile mixing are somewhat 
relaxed for high D/H. For low D/H (~ 2 x 10 -5 ), no active-sterile neutrino 
mixing is allowed by currently popular upper limits on the primordial 4 He 
abundance Y. For such low primordial D/H values, the observational infer- 
ence of active-sterile neutrino mixing by upcoming solar neutrino experiments 
would imply that Y has been systematically underestimated, unless there is 

new physics not included in standard BBN. 
14.60.Pq, 14.60.St, 26.35.+C 
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Upper limits on the abundance of 4 He produced in big bang nucleosynthesis (BBN) have 
been used to limit mixing between active (v e , u^, or z/ r ) and sterile (u s , no standard model 
interactions) neutrinos In this paper, we point out and discuss how these constraints 
are dependant on the adopted primordial deuterium abundance. Previous limits on sterile 
neutrino mixing have assumed a value for the lower bound on the baryon-to-photon ratio 
T) derived from interstellar medium and solar system measurements of deuterium (D) and 
3 He, and models of chemical and galactic evolution. Recent measurements of D/H in quasar 
absorption systems (QAS) have yielded discordant values of this ratio, some higher than 
previously derived ranges @, and some lower 0. Several factors make an investigation of 
the primordial D/H dependance of BBN constraints on active-sterile neutrino mixing timely: 
the discordant QAS measurements of D/H; the fact that future solar neutrino experiments 
may be able to distinguish and identify u e — v s mixing |4j]; and the use of sterile neutrinos 
in schemes for neutrino masses and mixings that explain all available data ||. 

As is well known (e.g., Ref. ||), the abundance of 4 He produced by BBN is essentially 
determined by the ratio of neutron to proton number densities (n/p) at "weak freeze-out" 
(WFO). WFO occurs when the reactions that interchange neutrons and protons proceed 
too slowly relative to the expansion rate of the universe to keep n/p at its equilibrium value 
of n/p « exp(— Am/T). Here Am = m n — m p m 1.293 MeV is the neutron-proton mass 
difference, and T is the photon temperature. Mixing between active and sterile neutrinos 
increases (n/p)wFO> an d therefore the primordial 4 He mass fraction Y, in two ways. First, 
active-sterile neutrino mixing effectively brings more degrees of freedom into thermal contact, 
increasing the energy density and hence the expansion rate of the universe. Second, active- 
sterile mixing — especially v e — v s mixing — depletes the electron neutrino and antineutrino 
populations, reducing the rates of the n <-» p interconversion reactions. Both of these effects 
cause n/p to freeze out at a lower temperature. 

Using a neutrino ensemble evolution formalism [I],[7| that includes both neutrino oscil- 
lations (with matter effects) and neutrino collisions, previous authors JIJ have produced 
exclusion plots in the Sm 2 -sm 2 29 plane for both v e — v s and — v s mixing ||. Here Sm 2 
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and sin 2 29 are the difference of the squares of the neutrino vacuum mass eigenvalues and 
a measure of the vacuum mixing angle, respectively, associated with two-flavor neutrino 
mixing. These studies showed that for 77 > 2.8 x 10~ 10 , both the — v s solution to the 



atmospheric neutrino problem [IIJ and the v e — v s large-angle MSW solution to the solar 
neutrino problem jD]] are excluded for Y < 0.247. 

In our study of the primordial D/H dependance of BBN constraints on active-sterile 
neutrino mixing, we have employed the same neutrino evolution formalism 0,^] as previous 
authors. We have neglected any net lepton number contributed by the neutrinos. (The 
recently reported effect of active-sterile neutrino mixing generating net lepton number does 
not occur in the regions of parameter space we consider here ||.) In this case, the neutrino 
and antineutrino sectors evolve identically. A Fermi-Dirac momentum distribution for all 
neutrinos is assumed, but allowance is made for non-equilibrium number densities. The 
differential equations in the formalism yield n Ue ,n u ^,n UT ,n UB , and T as functions of time. 
Here n Ux denotes the fraction of a full fermionic degree of freedom contributed by neutrino 
species x, which we shall hereafter call the "number density parameter" of neutrino species 
x. In the equations below we will take n Ux = n-p x , since we are working under the assumption 
that the net lepton number contributed by the neutrinos is negligible. 



In our BBN computation we have employed the Kawano |L2] update of the Wagoner |13 
code, with the latest world average neutron lifetime, r = 887.0 s fli~4| ; the reaction rates 
of Ref. jn|; and a correction of +0.0031 to Y due to finite nucleon mass and timestep- 



dependant effects fllEf . We have altered the Kawano code to use the 'temperature series' of 
neutrino number density parameters, n Ux (T), to compute the energy density contributed by 
neutrinos and the n <-> p interconversion rates. The neutrino energy density is 

7 7T 2 

Pv= ^-^(n UB + n Uti +n UT +n Us )T A . (1) 

The n <-» p rates are 

W = K f (j^) (l - l + n J { : +q)Zv ) x(x + g) V - 1) 1/2 dx, (2) 



W = K /~ ( l + n ;:_ q)Zv ) ( IT ^) x(x - g) V - l) 1 / 2 dz, (3) 

w = *f (y^i) (i - 1+e 7;^ ) *(* - <ov - !) 1/2 ^ 5 ( 4 ) 
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W = g/°°( 1 + " ( X^ ) (y^)^^^ 2 -!) 1 ' 2 ^, (6) 

W„ = A"/' (^) (j^^) x(» " «>V " l) 1/2 rf- (7) 

In these expressions a; = E e /m e , where E e and m e are the total electron (or positron) energy 
and rest mass, respectively; z = m e /T; z v = m e /T u , where is the appropriate neutrino 
temperature; q = Am/m e ; and X is a constant obtained by solving the equation 

for K, where r is the experimentally measured neutron lifetime. 

The lower limit on 77 obtained from a standard BBN calculation with N u = 3 is not 
appropriate for BBN with active-sterile mixing. This is because the lower limit on 77 depends 
on the expansion rate, often codified as an effective number of neutrino generations N u |17| . 
Since active-sterile mixing increases N v (at least for the range of parameter space we consider 
here ||), it affects the lower bound on rj. Therefore, we will plot our results as a function 
of the primordial D/H value — the experimentally determined quantity — rather than as a 
function of 77. These considerations are most important for the — v s atmospheric neutrino 
mixing solution, and much less important (nearly negligible) for the v e — v s small-angle MSW 
solution to the solar neutrino problem. 

In Fig. 1, a representative v^ — fs atmospheric neutrino mixing solution (5m 2 = 1.0 x 1CT 2 
eV 2 , sin 2 26 l = 0.6 |[J) is assumed, and the resulting BBN 4 He yield is plotted as a function 
of the BBN D/H yield. The value of 77 (given as 7710 = 10 10 rj) at various values of D/H is 



also indicated on the figure. For a given value of D/H, the implied abundance of 4 He can be 
interpreted as the observational upper limit required to constrain the solution. Alternatively, 
a detection of these neutrino mixing parameters by, for example, future atmospheric neutrino 
experiments would yield an independent determination of the primordial 4 He abundance, 
so long as D/H were known from QAS studies. This could be very interesting, given the 
recent emphasis on the systematic uncertainties in the determination of the primordial 4 He 
abundance as derived from helium recombination lines in extragalactic HII regions ||18|| . The 
conclusions reached from Fig. 1 are essentially the same over the range of 5m 2 (10~ 3 — 1CT 1 
eV 2 ) for the proposed — v s mixing explanation of the atmospheric neutrino problem. 

Fig. 2 is similar to Fig. 1, but with a representative v e — v s small-angle MSW solution 
to the solar neutrino problem assumed (5m 2 = 4.0 x 1CT 6 eV 2 , sin 2 26 l = 8.0 x 10~ 3 0). 
These mixing parameters have only a very small effect on the BBN 4 He yield. For the most 
recent data from solar neutrino experiments and the standard solar model p0| , there is no 
v e — v s large-angle MSW solution to the solar neutrino problem [|HJ . 



Some of the QAS data suggest D/H « 2 x 10~ 4 0. Figs. 1-2 show the range D/H = 
1.5— 2. 3x 10~ 4 , as determined in Ref. [21]. This range of D/H implies a lower bound on rj that 



is significantly lower than that used in previous studies. Since a lower 77 implies a lower 4 He 
yield, high D/H relaxes constraints on any effect that increases the expansion rate, including 



mixing with sterile species [22] . Fig. 1 shows, however, that for D/H « 2 x 10 , the v„ — v s 



atmospheric neutrino solution is still somewhat constrained [^3] if current observational 
inferences [p4 of primordial 4 He are correct: Y = 0.234 ± 0.003 ± 0.005, where the first 



error is statistical and the second systematic. Of course, if this z/ M — v s atmospheric neutrino 
mixing solution were inferred from atmospheric neutrino experiments, and QAS studies 
confirm D/H « 2 x 10~ 4 , the implied 4 He abundance of Y w 0.245 would be significantly 
higher than the central value of Y = 0.234 cited above. 

The v e — v s small-angle MSW solar neutrino solution is allowed for D/H w 2 x 10 -4 . Fig. 
2 shows that an observational upper bound of Y ^ 0.232 would be required to restrict this 
small angle solution if such a high D/H is indeed the primordial value. 



Other very high quality QAS data — arguably better |3j for the determination of D/H 
than that used in Ref. H — suggest D/H k2x 10~ 5 0]. This value of D/H is incompatible 



with standard BBN with N v = 3 |Tq|25|fi7H for current observational inferences of the primor- 
dial 4 He abundances |24|], and any mixing with sterile neutrinos would only exacerbate the 
problem. As mentioned previously, however, it has been argued that Y has been systemati- 
cally underestimated, and a more appropriate upper limit on Y may actually be Y < 0.255 
|18| . It is unlikely that the systematic error in Y is enough to allow the — v s atmospheric 



neutrino solution for D/H rs 2 x 10~ 5 . However, observation of the u e — v s small-angle MSW 
solar neutrino solution, together with a solid determination of D /H « 2 x 10~ 5 , would require 
that Y has been systematically and significantly underestimated by about 0.015 (see Fig. 2), 
unless there is non-standard physics during the BBN epoch [p5| |. This is a somewhat trivial 



point, since the mixing parameters of the small-angle v e — v s MSW solution produce only 
very slightly more 4 He than the standard BBN picture with N u = 3, for which the "crisis" 



at low D/H is well-known pl||2"5| , |T7f . Useful constraints on the u e — v s small-angle MSW 
solar neutrino solution would require very precise observational knowledge of r] and Y. This 
may, however, still be interesting in view of the fact that future solar neutrino experiments 
may be able to distinguish the sterile neutrino oscillation-based solution from other solutions 
Also, many models that seek to satisfy all available constraints on neutrino properties 



employ the u e — v s small-angle MSW solar neutrino solution || (but see Ref. [26] ) 
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FIGURES 

FIG. 1. BBN yields for a typical — v s atmospheric neutrino solution (5m? = 1.0 x 10~ 2 eV 2 , 
sin 2 26 = 0.6). The solid curve is the 4 He mass fraction Y vs. D/H. The squares indicate, from 
lower left to upper right, 10 10 r] = 1.7, 2.3, 3.0, 4.6, 6.6, 8.6. The dotted lines indicate the ranges 
of "high" and "low" D/H inferred from QAS studies. The dashed lines indicate two possible upper 
limits on Y: the currently popular Y = 0.245, and the more conservative Y = 0.255. 

FIG. 2. BBN yields for a typical v e — v s small-angle MSW solution to the solar neutrino problem 
(,5m 2 = 4.0 x 10~ 6 eV 2 , sin 2 26 = 8.0 x 10~ 3 ). The solid curve is the 4 He mass fraction Y vs. D/H. 
The squares indicate, from lower left to upper right, 10 10 -q = 1.5, 2.0, 2.6, 3.6, 4.6, 5.6. The 
dotted lines indicate the ranges of "high" and "low" D/H inferred from QAS studies. The dashed 
lines indicate two possible upper limits on Y: the currently popular Y = 0.245, and the more 
conservative Y = 0.255. 
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